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Full-potential gradient corrected density functional calcu- 
lations of the adsorption of potassium on the Ag(lll) sur- 
face have been performed. The considered structures are 
Ag(lll)(V3x V3)R30°-K and Ag(lll)(2x2)-K. For the lower 
coverage, fee, hep and bridge site; and for the higher coverage 
all considered sites are practically degenerate. Substrate rum- 
pling is most important for the top adsorption site. The bond 
length is found to be nearly identical for the two coverages, 
in agreement with recent experiments. Results from MuUiken 
populations, bond lengths, core level shifts and work func- 
tions consistently indicate a small charge transfer from the 
potassium atom to the substrate, which is slightly larger for 
the lower coverage. 



I. INTRODUCTION 

The study of adsorbates on metal surfaces has become 
a center of interest in surface science because of the enor- 
mous importance of catalysis for industrial applications. 
These experimental studies can be complemented with 
simulations which are nowadays a powerful tool for the 
theoretical description of surface structures and chemical 
reactions at surfaces, because of the growth in computa- 
tional power and due to the success of density functional 
calculations. 

One prototype reaction is the adsorption of alkali met- 
als on metallic surfaces. Although this is a relatively 
simple reaction, there have been several surprises in the 
study of these systems. Initially, it was assumed that the 
adatoms would occupy hollow sites, until a first system 
was discovered where top site adsorption was preferredEI. 
Meanwhile, further systems of alkali metals with top site 
adsorption have been found — typically heavier alkali 
atoms (from K on) on close packed surfaces (for reviews, 
see, e.g., references |[ |[ |[ ^, 0) . One important find- 
ing was that top site adsorption was accompanied by 
substrate rumpling and thus it was assumed that sub- 
strates with a low bulk modulus would be prominent 
for top site adsorption. It was therefore a little surprise 
when in a LEED (low energy electron diffraction) study 
of K on Ag( 111) it was discovered that hollow sites were 
occupiecotl. 

The type of hollow site occupied in the system 
K/Ag(lll) depends on the coverage: in the low cov- 
erage structures ((2 x 2) and (3x3)), the adatoms 
occupy face centered cubic (fee) hollow sites, in the 
(-s/S X a/3)R30° structure the hexagonal close packed 



(hep) hollow sitefii. In addition, the phase diagram 
of K, Rb and Cs adsorbed on Ag(lll) was investigated 
exp er iment allyll3 . 

There have been several simulations of alkali metals on 
close-packed metallic surfaces: for example, systems such 
as Na fuid K on the Al(lll) surfaceliil, Na on the Cu(llli 
surfaced, K on the Pt(lll) surfacell^and K on Cu(lll)Eil 
have been studied. The system K/Ag(lll) is therefore a 
very interesting system due to the experimental findings 
and as an extension of the employed technology (section 
H) to substrates of the second row of the transition met- 
als. 

It is therefore the aim of this article to summa- 
rize results from simulations on the system K/Ag(lll). 
Full-potential density-functional calculations with a local 
Gaussian basis set and with a gradient corrected func- 
tional were performed, for two coverages (1/3 or 1/4 of a 
monolayer, for a (V3x V3)R30° or (2 x 2) pattern). The 
addressed questions are the preferred adsorption sites, 
the magnitude of the energy splitting of different highly 
symmetric adsorption sites, the geometry and the impor- 
tance of surface rumpling, the charge of the adsorbate, 
the work function and the positions of the K core levels. 



II. COMPUTATIONAL PARAMETERS 

A local basis set formalism was used where the ba- 
sis functions are Gaussian type orbitals centered at the 
atoms as implemented in the code CRYSTALEj. For Ag, 
a relativistic small-core pseudopotential (with 28 elec- 
trons in the core) in combination with a [4s3p2(i] basis set 
was employed as in a previous study of CI on Ag(lll)E3. 
The [5s4pl(i] K all-electron basis from reference |l4| was 
chosen. The exchange-correlation potential fHf^ fitted 
with auxiliary basis sets as in previous worktil'tj. The 
gradient corrected exchange and correlation functional 
of Perdew, Burke and Ernzerhof (PBE)0 was used. 

The adsorption was modeled by using slabs of 4 or 
5 silver layers, at the PBE-optimized bulk Ag lattice 
constantllj of ao = 4.10 A. Potassium was adsorbed 
on one side of this slab. A supercell approach with a 
(V3 X V3)R30° or (2 x2) structure as in the experiment 
was implemented. This slab was truly two dimensional 
and thus not periodically repeated in the third dimen- 
sion. The vertical relaxation of the silver atoms was 
simulated in three different ways: simulations were per- 
formed where only a uniform relaxation of the top silver 
layer was possible, and simulations where a different ver- 
tical relaxation of the silver atoms in the top layer was 
possible (i.e. substrate rumpling). Finally, an additional 
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lateral displacement of the atoms in the tcpiSilver layer, 
as observed experimentally for K/Ni(lll)E3, was simu- 
lated. 

Four adsorption sites were considered (see figures ^ and 
H) : the top adsorption site with K sitting vertically above 
a silver atom in the top layer, the bridge site with K 
sitting above the middle of two silver atoms in the top 
layer, and two different threefold hollow sites where the 
potassium atoms are placed vertically above a silver atom 
in the second (third) silver layer (hep and fee hollow, 
respectively). In addition, the possibility of a stacking 
fault was investigated. In this structure, the outermost 
Ag layer and the K layer would occupy hep sites. Such 
a stacking fault was observedll3 for the system Ag(lll) 
(73 X \/3)R30°-Sb. 

The optimized geometrical parameters are defined as 
in figure ^. The variables are chosen in such a way so 
that the notation igjConsistent with those from the exper- 
imental referencesErLI. The bond length is defined as the 
distance from a potassium adatom to the nearest silver 
atom. dx-Agi, dAgi-Ag2, dAg2-Ag3 and SI are various 
interlayer distances: dK-Agi is the vertical distance from 
the potassium layer to the plane made of the atoms in the 
first silver layer which have moved towards the potassium 
layer. dAgi-Ag2 is the distance between the plane made 
of the atoms in the first layer which have moved towards 
the second Ag layer, and the second Ag layer. dAg2-Ag3, 
the distance between second and third layer silver atoms, 
was kept fixed at the value of 4.10/-\/3 A corresponding 
to the distance of two layers in the bulk (ao / V^) ■ Finally, 
(51 is the rumpling, i.e. the distance between the plane 
made of those atoms in the first silver layer which have 
moved towards the K atoms and the plane made of those 
atoms in the first silver layer which have moved towards 
the second silver layer. Additionally, the possibility of a 
lateral displacement of the atoms in the top Ag layer was 
taken into account (for a more detailed description of the 



lateral displacement, see the following section [II 



A net with 16 x l|6-Bampling points in the surface Bril- 
louin zone was usedE3. During the optimization process, 
the Fermi function was smeared with a temperature of 
0.03 Eh (l£',i=27.2114 eVX and the total energy extrapo- 
lated to zero temperature^. To ensure the stability of the 
results with respect to this parameter, additional single- 
point calculations with a smearing temperature of 0.005 
Eh were performed (see also appendix). Also, the results 
for properties such as MuUiken population, work function 
etc. were computed with a smearing of 0.005 Eh- For an 
extensive test of the various computational parameters 
for metallic systems with the CRYSTAL code, see also 
reference 

The unit cells which are used in the simulations for 
this article are at the limit of what is presently feasible 
with this method. This local basis set approach is thus 
for metals probably more expensive than alternative ap- 
proaches relying on plane wave technologies. However, it 
is interesting to have the possibility of comparing various 



strategies which complement each other. One advantage 
of the method employed here is that concepts such as 
computing MuUiken charges are easily implemented in a 
code relying on Gaussian basis sets. Also, all-electron 
calculations do not pose a problem, and although in the 
present work a small-core pseudopotential was used for 
Ag, an all-electron description would also have been pos- 
sible. The usage of a pseudopotential has the advantage 
that scalar-relativistic corrections can be included. 



III. RESULTS AND DISCUSSION 

A. Ag(lll) {VS X V3)R30°-K 

In this section, the results of the calculations on the 
system Ag(lll) {V^x •\/3)R30°-K are presented and dis- 
cussed. 

In table |, the adsorption geometries and energies are 
displayed. It turns out that all considered adsorption 
sites are virtually degenerate — the best calculations in- 
dicate that the binding energies for the various sites are 
within 0.1 mEh which is at the limit of the numerical 
noise. These findings are thus no contradiction to the 
experiment where the hep hollow site was identified as 
the preferential site. The degeneracy is certainly within 
the range of the errors associated within the calculation 
(lack of basis functions with higher angular momentum 
than 1 = 2, number of Ag layers considered in the slab 
model, uncertainties in the functional, anisotropic vibra- 
tions observed in the experiment). The total computed 
binding energy per K atom is 0.042 Eh, i.e. 1.1 eV. 

The structure with a stacking fault is also degenerate, 
and thus cannot be ruled out bypthe simulations. In 
agreement with a recent calculatiorS, the clean Ag(lll) 
surface with a stacking fault was found to be slightly 
higher in energy (by about 0.00015 Eh per atom) than 
the unfaulted surface. 

When considering the optimal geometry for the hep 
site, we note that the computed distances are in excel- 
lent agreement with the experiment: the distance be- 
tween the potassium layer and the top silver layer dx-Agi 
was found rto be 2.79 A versus 2.84 ± 0.03 A in the 
experiments and the distance between first and second 
silver layer dAgi-Ag2 was obtained to be 2.37 A (exp.: 
2.35 ± 0.02 A). This computed distance is identical to 
the distaBee in Ag bulk (4.10/\/3 A, with the PBE 
functionalij) , i.e. no relaxation of the top layer was 
found. The distance between potassium and the near- 
est silver atom was found to be 3.27 A (exp.: 3.29 ± 0.02 
A). We can thus deduce an effective potassium radius of 
3.27— ^^A=1.82 A, in agreement with the experimental 

valuJ of 1.82 ± 0.03 A. 

The most complex structure included the possibility of 
surface rumpling and of lateral relaxation of the top Ag 
layer. Rumpling is allowed by symmetry for the bridge 
and top adsorption site, and in both cases a large value 
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51 for the rumpling is obtained (0.10 A for the bridge site 
and 0.22 A for the top site). The effect of the rumpling 
is to push the Ag atom under the adatom deeper into the 
bulk, whereas the atom(s) not under the K atom move 
out of the bulk towards the K layer. 

The importance of rumpling in the second Ag layer un- 
der the K layer was investigated for the hep hollow site, 
where this is allowed by symmetry. An optimal displace- 
ment of 0.01 A of the Ag atom vertically under the K 
atom, away from the K atom, was found, with an energy 
gain of 0.1 mE^. In the experiment, a displacement of 
0.02 ± 0.02 A towards the K adsorbate layer was found 
for the Ag atom in the second layer under the K atom. 
Although the computed result deviates thus slightly from 
experiment, it is important that the energy gain is fairly 
small so that second layer rumpling does not appear to 
have a huge influence on the results (this was suggested 
in reference [s] as a possible reason for the switch in ad- 
sorption site). 

Lateral displacement in the top Ag layer plays only a 
minor role, and lowers the energy at most by 0.4 mEfi. 
The atoms nearest to the K adsorbate were allowed to 
move and the lateral movement of the atoms was in all 
cases away from the adsorbate, by a maximum of 0.04 A. 

We also note by comparing the results for 4 layers with 
and without surface rumpling, that the rumpling is most 
important for the top site: the energy gain by rumpling 
is 2.5 mEh, and 0.5 mEh for the bridge site. .-This is in 
line with previous calculations for K/Cu(lll)EJ. 



B. Ag(lll) (2 X 2)-K 

In the following table ||, results for the Ag(lll)(2 x 2)- 
K structure are displayed. Fee hollow, hep hollow and 
bridge site are within 0.2 mEh, i.e. practically degener- 
ate. The structure with a stacking fault is ^ 0.8 mEh 
higher in energy, the top site by ~ 1.8 mEh- Similar as 
for the higher coverage, it can only be stated that these 
findings are compatible with the experiment, but a state- 
ment about the preferred site is not possible. 

The geometry is in good agreement with experiment: a 
value for the interlayer distance dx-Agi of 2.60 A (exp.: 
2.70 ± 0.03 A) was obtained and an interlayer distance 
dAgi-Ag2 between the first and second silver layer of 2.35 
A (exp.: 2.34 ± 0.02 A). Surface rumpling is important 
and a value of 61 =0.11 A was computed (exp.: 0.10 
± 0.03 A). The effect of the surface rumpling was again 
to push those atom(s) under the potassium adsorbate 
deeper into the substrate, whereas the other atom (three 
atoms in the case of the top site, two atoms in the case of 
the bridge site and one atom in the case of the threefold 
hollow sites) move outwards towards the potassium layer. 

Lateral relaxation plays a minor role for the binding 
energy. In all cases except for the top site where this 
is not possible, the atoms nearest to K move away from 
the K atom; in the case of the top site, the three atoms 



not under the K atom move slightly towards each other, 
away from the K atom. However, the energy gain is at 
the limit of the numerical accuracy. 

For comparison, calculations without surface rumpling 
were performed, i.e. 61 was kept fixed at A. We note 
that the interlayer distances dx-Agi between K layer and 
top Ag layer increase when rumpling is not possible, and 
are 2.73 A for fee and hep site. The interlayer distance 
between first and second Ag layer hardly changes. The 
adsorption energy differs by < ImEh for fee, hep, bridge 
site and for the structure with a stacking fault; but by 2.9 
mEh for the top site. Thus we find again that surface 
rumpling is most important for the top site, however, 
it is not big enough to make the top site the preferred 
adsorption site. 

The computed adsorption energy per potassium atom 
for this structure is 0.041 Eh, i.e. nearly identical to the 
higher coverage. This energy is close to the one com- 
puted for K/A1(111); also the splitting between various 
sites was fojujid to be of the order of few hundreths of 
an eV onlyEiJ; thcj-binding energy is also comparable to 
Cu(lll)(2 X 2)-KEj. The adsorption energy is however 
lower than the one computed for K/Pt(lll) where values 
of 2.42 eV (at a coverage of one third of a monolayer) 
and 2.93 eV (at a coverage of one fourth of a mono- 
layer) were obtained, at the level of the local density 
approximationlj. The energy splitting between fee and 
hep site was computed to be 15 meV in the latter publi- 
cation, with the hep site being lower; this is compatible 
with the near-degeneracy which is found here for the sys- 
tem K/Ag(lll). 

Finally, we also see that the results are well converged 
with respect to the number of layers, for both coverages: 
in one case (substrate rumpling, no lateral relaxation), 
the optimization has been performed with 4 and 5 Ag 
layers, and both geometry and adsorption energies are 
stable. 



C. Population and density of states 



In table III, MuUikcn charges for the potassium atom 



are given, projected on the different basis functions. 
These charges are defined as 



with density matrix P, overlap matrix S and A the set 
of basis functions for which the population is computed. 

It should be mentioned that a relatively large number 
of digits has been given; although the absolute values will 
not be more accurate than to one or two decimal points, 
for the relative values more digits are still meaningful. 

In the case of Ag(lll) (V3 x V3)R30°-K, the total 
charge is ^ 0.16 |e|, in the case of Ag(lll) (2 x 2)-K 
^ 0.24 |e|. Obviously, if the charges were the same, 
the electrostatic repulsion at higher coverage (Ag(lll) 
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(-\/3 X -\/3)R30°-K) would be larger because of the shorter 
K-K distances. Therefore, the charge in the case of 
Ag(lll) (-s/Sx -\/3)R30°-K is reduced, in agreement with 
the Langmuir-Gurney model. The individual populations 
are very similar to Cu(lll) (2 x 2)-K: p^, Py and Pz 
charges are slightly above 4, the number of s-electrons 
is 6.5; charges in p^ and Py orbitals are identical ex- 
cept for the bridge site; the Px orbital which has more 
overlap with the substrate layer (in our choice of geom- 
etry) is slightly less occupied than the Py orbital. The 
higher occupancy of that orbital, which has less overlap 
with the substrate, indicates that the bond is not cova- 
lent. In addition, the overlap population defined as 

with A the set of basis functions on the first atom 
and B the set of basis functions of the other atom, was 
computed. This number is a measure for covalency, and 
is, similar to Cu(lll)(2x2)-K, very small: for K and 
nearest neighbor Ag, it is < 0.06 |e| for both coverages, 
and negligible for further neighbors. 

In figures ^ and ||, the density of states (DOS), pro- 
jected on the potassium basis functions, is displayed for 
the fee adsorption site. The Fermi energy is in a regime 
with a significant contribution from K 4s and 4p bands, 
so that the overlayer is clearly metallic. The projected 
DOS does not depend on the adsorption site: it looks 
virtually identical for fee hollow, hep hollow, bridge and 
top site, and also for the structure with stacking fault. 
This is as found for Cu(lll)(2 x 2)-K, but different from 
the case of chlorine as an adsorbate, where the-projected 
DOS clearly depended on the adsorption sitell3. 

The values of the K 3s and 3j5 levels are practically 
independent of the adsorption site (see table IV). This 
is again in contrast to chlorine: for chlorine, the core 
eigenvalue depended on the adsorption site, and-it was 
correlated with the charge of the chlorine atomlla. For 
K, the variation of the MuUiken charges with adsorption 
site for fixed coverage is much smaller (table III) as for 
chlorine. This is thus consistent with the finding that the 
core eigenvalues do not vary for the different sites. Also, 
the position of the Fermi energy and thus the value of 
the work function is independent of the site, and changes 
slightly with the coverage. 



D. Comparison of the results for the two coverages 

When comparing the results for the two coverages, we 
firstly note that the bond lengths for the individual sites 
are very similar for both coverages. The bond length in- 
creases with the number of nearest neighbors to which 
the potassium jadatom bonds, in agreement with Paul- 
ing's argumentcj that a smaller number of bonds leads 
to a stronger individual bond and a shorter bond length: 
the bond length is shortest for the top site with 2.95 A 
for the (V3 x %/3) structure (or 2.89 A for the (2 x 2) 
structure), longer for the bridge site and longest for the 



threefold hollow sites. Similarly, the effective K radius, 
obtained from the difference of bond length and effective 
silver radius (^^A), increases in this order. If we keep 

the Ag radius fixed, it increases from 1.50 A (or 1.44 A 
for the (2 x 2) structure) for the top site to 1.82 A (or 1.76 
A) for the threefold hollow sites. These values are com- 
patible with effective K radii deduced from experiment 
(see table 10 in reference 

In table 0, the most important results of the calcula- 
tions for the two coverages and the fee hollow arc com- 
pared. 

The binding energy per K atom, with respect to free 
K atoms and a clean Ag(lll) slab, is 0.042 and 0.041 
Eh, i.e. nearly identical for both coverages. In the case 
of Ag (111) (^/3 X V3)R30°-K, aU considered sites, and 
in the case of Ag(lll) (2 x 2)-K, the fee, hep and bridge 
site are virtually degenerate. This degeneracy was al- 
ready observed in the early study of K/Al(lll)Eil and 
explained with the fact that the large radius of potas- 
sium leads to a large distance of the adsorbate to the 
substrate and therefore the K adsorbate will experience 
only a small substrate electron density corrugation; simi- 
larly the potassium overlayer is metallic and its electmnic. 
charge mpce diffuse than the charge of, e.g. chlorideg£3'E£l 
or oxidesEj as adsorbates which also helps to reduce the 
energy difference between the threefold hollow sites and 
bridge or top site, compared to halogenides or oxides. 

The K charge hardly depends on the adsorption site, 
but it depends on the coverage. The slightly higher 
charge (i.e. less electrons) for the lower coverage is con- 
sistent with the slightly shorter bond length and smaller 
effective K radius for this adsorbate system. 

The values of the K 3s and ip core eigenvalues also 
indicate that there is a little difference between the two 
coverages: For the lower coverage, the core eigenvalues 
are at 1.194 Eh and 0.590 Eh below the Fermi level, for 
the higher coverage at 1.184 Eh and 0.600 Eh- These 
values are virtually independent of the adsorption site. 
Again, this is consistent with the finding that the K 
charge is larger for the lower coverage and therefore the 
core eigenvalues are lower (i.e. there is slightly less elec- 
tronic charge to screen the nuclear charge and the core 
levels are thus further stabilized). 

The work function decreases from 0.131 Eh for the 
clean Ag(lll) surface to 0.062 Eh for the (2 x 2) struc- 
ture and increases again slightly to 0.069 Eh for the 
(V3 X V3)i?30 structure. This finding is in line with the 
argument that initially, the work function will decrease 
linearly with increasing coverage, up to a minimum, and 
finally increase again because of the depolarization of the 
adsorbate (see, e.g. ref. ||). 

Comparing the results for K/Ag(lll) with K/Cu(lll) 
and the (2 x 2) pattern, we note that the site has 
switched from top to the threefold hollow or bridge site. 
The reason appears to be that, without rumpling, the 
energy splitting between hollow or bridge sites and top 
site is less than 1 mEh for K/Cu(lll)Ej, but 4 mEh for 
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K/Ag(lll). Substrate rumpling changes this splitting in 
favor of the top site; however, the change is not large 
enough to stabilize the top site as the ground state for 
K/Ag(lll). The parameter most important for the mag- 
nitude of the energy splitting of the sites, without rum- 
pling, appears to be the nearest neighbor spacing of the 
substrate: for a fixed radius of the K adsorbate, the ad- 
sorbed atom will have more overlap with more substrate 
atoms, if the substrate atoms are closer. In the case of 
the top site, the K atom has a nearest neighbor vertically 
below and six next nearest neighbors in the first layer. 
These six atoms are closer in the case of K/Cu(lll) than 
K/Ag(lll) because of the smaller atomic radius of cop- 
per and thus smaller lattice constant of Cu. Thus, the 
energy splitting is smaller for K/Cu(lll) compared to 
K/Ag(lll), and substrate rumpling can change it in favor 
of the top site for K/Cu(lll). Tius argument would also 
hold for the system K/Ru(0001)E3 which behaves similar 
to K/Ag(lll). The calculations thus support the argu- 
ment given, for example, in reference |^ and which was 
based on a large set of experimental data (when varying 
the adsorbate, the ratio of the radius of the adsorbate to 
the nearest neighbor spacing of the substrate was recom- 
mended as a parameter to predict the adsorption site). 
In the simulations performed here, this argument can be 
verified by comparing the results from a surface without 
substrate rumpling and a surface with rumpling. 

For substrates with similar lattice constants, an addi- 
tional parameter will be the energy required for the defor- 
mation of the surface; this energy is found to be 1.5 mEii 
for K/Cu(lll) and 2.8 mEh for K/Ag(lll), for the top 
site ((2 X 2) structure, the energy difference between the 
unrelaxed clean surface and the optimal geometry of the 
adsorbate system, without K adsorbate, was computed). 
Note that the magnitude of the rumpling is larger for 
K/Ag(lll), which may explain the higher energy neces- 
sary for the deformation of the clean Ag surface (although 
Ag has a lower bulk modulus than Cu). 

We also note that the top site is slightly better in en- 
ergy for the higher coverage than for the lower cover- 
age, relati-pi to the other sites (in line with the case of 
K/Al(lll)y). This could be due to the larger radius 
of the K adsorbate for the higher coverage. As explained 
earlier, the positive K charge is smaller for the higher cov- 
erage and thus the radius is larger. Therefore, the overlap 
with the 6 nearest neighbors will be larger and the top 
site will be slightly more stabilized. However, this argtu 
ment must be contrasted with the case of Cs/Ru(0001)lZI 
where the site switches from top ((2 x 2) structure) to 
hep ((%/3 X %/3)R30° structure). In the latter case, a 
better screening was suggested as the possible explana- 
tion. These two controversial findings for similar systems 
demonstrate how difficult it is to give even only qualita- 
tive rules for the preferred adsorption site. 

We can thus give a reasonable argument for the ques- 
tion, why the top site is not preferred for K/Ag(lll), in 
contrast to K/Cu(lll). The question of the site switch 



from fee to hep, however, can not yet be answered. The 
difficulty results from the small energy splitting, and thus 
large number of effects on this energy scale that are pos- 
sible, and may influence the site preference (e.g. slightly 
different charge, screening, second layer rumpling, stack- 
ing fault). Numerically, the energy splitting is at the 
limit of the numerical noise. 



IV. SUMMARY 

In summary, the impact of the article is twofold: firstly, 
the technology of using Gaussian type basis functions 
was shown to be capable of describing the adsorption 
on metallic surfaces. Secondly, some of the experimen- 
tal findings could be confirmed (essentially the geome- 
try), and new data could be provided (energy splitting 
between the various adsorption sites, Mulliken charges, 
core eigenvalues, work functions, and the dependence of 
these parameters on the coverage). 

In more detail, it has been shown that density func- 
tional calculations can reproduce a part of the experi- 
mental findings for the systems Ag (111) (a/S x •\/3)R30°- 
K and Ag(lll) (2x2)-K. The computed geometry is in 
excellent agreement with the experimental results, and 
therefore much better than for the case of K/Cu(lll) 
where there was a deviation QjLpi.0.2 A between the ex- 
periment and the simulationEao. This agreement for 
K/Ag(lll) thus helps to support the validity of the ap- 
proach employed in the simulations. 

The adsorption site was found to be nearly degener- 
ate for various sites considered for both coverages, and a 
value of 0.04 Eh was obtained for the binding energy. The 
main difference to the system K/Cu(lll), where the top 
site is occupied, appears to be due the different lattice 
constant of Cu vs. Ag. In both threefold hollow 

site would be favored without substrate rumpling, but 
only in the case of K/Cu(lll), the energy gain associated 
with the rumpling is large enough to lead to a change of 
the adsorption site and makes the top site favorable. 

The computed energies demonstrate that various sites 
are virtually degenerate, within the accuracy of the cal- 
culations. This appears to be due to the metallic nature 
of the overlayer and the large radius of the potassium 
adsorbate, which makes the surface electron-density cor- 
rugation have only little impact on the energy splitting 
between the various sites. Possible suggestions for the 
change in adsorption site from fee to hep were tested 
(second layer rumpling, stacking fault). However, sec- 
ond layer rumpling did not have a huge impact, and the 
energy of the structure with a stacking fault was still de- 
generate with the other structures, within the accuracy 
of the calculations. It is therefore possible to give argu- 
ments why the top site is not occupied, but the question 
of the site change from fee to hep remains puzzling. The 
near-degeneracy of the energies indicates that minute ef- 
fects are responsible for the site change. 
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The computed potassium charge is between +0.16 |e| 
for the higher and +0.24 |e| for the lower coverage. These 
charges are consistent with results for the effective K ra- 
dius and K-Ag bond length which increase for the higher 
coverage, with the work function, and with the K 3s 
and 3p core eigenvalues which move slightly closer to the 
Fermi energy for the higher coverage. This consistency 
indicates that MuUiken charges are a simple and useful 
concept for describing adsorption on metallic surfaces. 
As in the case of K/Cu(lll), the small overlap popula- 
tion indicates that virtually no evidence can be found for 

a covalent contribution to the binding. 

* 



APPENDIX A: NUMERICAL INTEGRATION 

The accuracy of the numerical integration depends es- 
pecially on the number of sampling points and the smear- 
ing temperature employed. Obviously, a higher number 
of sampling points will lead to a higher accuracy. The 
optimal smearing temperature will depend on the band 
structure and on the number of fc-points. In general, a 
lower smearing temperature will be better because the 
Fermi function (which is a step function at zero temper- 
ature and thus makes accurate integration numerically 
difficult) is better approximated. However, when quan- 
tities such as susceptibilities are computed, the smearing 
temperature should not be too low, because the function 
to be integrated has singularities (see, e.g., the discussion 
in reference pT] ). 

To investigate the dependence of the results on the 
smearing temperature, first a simple test was performed 
where the lattice constant of bulk Ag was optimized at 
T = O.OObEh- The optimal lattice constant was 4.12 A, 
the cohesive energy 0.090 Eh and the bulk modulus 108 
GPa. At a higher J^mperature of T = 0.03Eh, the opti- 
mized values wereEEl 4.10 A, 0.088 Eh and 113 GPa. The 
higher smearing temperature of 0.03 Eh thus does not 
seem to have a huge impact on the result. This is also 
seen in the comparison in tables | and where the en- 
ergy splitting between the various sites hardly changes 
when the smearing temperature is reduced from 0.03 Eh 
to 0.005 Eh- There is a change in the total binding en- 
ergy, which is however still negligible compared to its 
magnitude. 

Finally, in figure |, the binding energy Ex/Ag{iii) ~ 
Ek — EAg(iii) is computed for the hep site, with a 3 
layer slab simulating the Ag(lll) surface. It becomes 
obvious that the binding energy depends similarly on the 
number of fc-points as on the smearing temperature. As 
a whole, it can be concluded that the employed smear- 
ing temperature should not lead to serious errors for the 
results presented in this article. 
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TABLE I. Adsorption of K on the Ag(lll) surface, {\/3 x ^/3) R30° pattern. The adsorption energy Eads is the difference 
^K/Ag(iii) — ^'Ag{iii) — Ek- For the most complex geometry, this energy was computed with two different parameters of the 
smearing temperature, 0.03 Eh and 0.005 Eh, for comparison. 

site Bond length dx-Agi dAgi-Ag2 SI lateral Eads{0.03) -Eads (0.005) 





[A] 


[A] 


[A] 


[A] 


[A] 


r E, 1 

L J\ dto TYl- \ 


r E,, ■ 

L CL tOTfZ 






4 layers, without substrate rumpling, without lateral relaxation 






fee hollow 


3.26 


2.80 


2.37 








-0.0442 




hep hollow 


3.26 


2.80 


2.37 








-0.0440 




bridge 


3.16 


2.81 


2.37 








-0.0437 




top 


2.94 


2.94 


2.37 








-0.0413 




stacking fault 


3.27 


2.81 


2.38 








-0.0435 








4 layers, with substrate rumpling, without lateral relaxation 






fee hollow 


3.26 


2.80 


2.37 








-0.0442 




hep hollow 


3.26 


2.80 


2.37 








-0.0440 




bridge 


3.18 


2.73 


2.34 


0.10 





-0.0442 




top 


2.96 


2.73 


2.24 


0.23 





-0.0438 




stacking fault 


3.27 


2.81 


2.38 








-0.0435 








5 Ag layers, with substrate rumpling, without lateral relaxation 






fee hollow 


3.27 


2.81 


2.37 








-0.0443 




hep hollow 


3.27 


2.81 


2.37 








-0.0441 




bridge 


3.18 


2.74 


2.34 


0.10 





-0.0443 




top 


2.95 


2.73 


2.25 


0.22 





-0.0441 




stacking fault 


3.27 


2.81 


2.38 








-0.0437 








5 Ag layers, with substrate 


rumpling 


with lateral relaxation 






fee hollow 


3.27 


2.79 


2.37 





0.04 


-0.0444 


-0.0422 


hep hollow 


3.27 


2.79 


2.37 





0.04 


-0.0442 


-0.0421 


bridge 


3.19 


2.74 


2.33 


0.10 


0.01 


-0.0444 


-0.0422 


top 


2.95 


2.73 


2.25 


0.22 





-0.0441 


-0.0422 


stacking fault 


3.27 


2.79 


2.38 





0.03 


-0.0441 


-0.0422 



exp (hcp)il 3.29 ± 0.02 2.84 ± 0.03 2.35 ± 0.02 
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TABLE II. Adsorption of K on the Ag(lll) surface, (2 x 2) pattern. The adsorption energy Eads is the difference 
i?K/Ag(iii) ^ i-AgCiii) — Ek- For the most complex geometry, this energy was computed with two different parameters of 
the smearing temperature, 0.03 Eh and 0.005 Eh, for comparison. 



site Bond length dK-Agi dAgi-Ag2 51 lateral Eads{O.OS) Eads{0.005) 





[A] 


[A] 


[A] 


[A] 


[A] 


r 1 

L K atom J 


r E,, ■ 

I K atom . 






4 Ag layers, without substrate rumpling 


I, without lateral relaxation 






fee hollow 


3.20 


2.73 


2.37 








-0.0418 




hep hollow 


3.20 


2.73 


2.37 








-0.0416 




bridge 


3.11 


2.75 


2.37 








-0.0412 








2.87 


2.36 








U.UO ( o 




stacking fault 


3.21 


2.74 


2.39 








-0.0409 








4 Ag layers, with 


substrate rumpling, without lateral 


relaxation 






fee hollow 


3.21 


2.61 


2.34 


0.13 





-0.0424 




hep hollow 


3.21 


2.61 


2.34 


0.13 





-0.0422 




bridge 


3.12 


2.64 


2.31 


0.12 





-0.0422 




top 


9 on 


2.67 


2.22 


0.23 









stacking fault 


3.21 


2.61 


2.35 


0.13 





-0.0415 








5 Ag layers, with 


substrate rumpling, without lateral 


relaxation 






fee hollow 


3.20 


2.61 


2.34 


0.12 





-0.0423 




hep hollow 


3.21 


2.61 


2.34 


0.13 





-0.0423 




bridge 


3.12 


2.66 


2.32 


0.11 





-0.0422 




top 


2.89 


2.66 


2.22 


0.23 





-0.0408 




stacking fault 


3.21 


2.62 


2.35 


0.12 





-0.0415 








5 Ag layers, with substrate rumpling 


with lateral relaxation 






fee hollow 


3.20 


2.60 


2.35 


0.11 


0.04 


-0.0424 


-0.0411 


hep hollow 


3.21 


2.60 


2.35 


0.12 


0.03 


-0.0425 


-0.0413 


bridge 


3.13 


2.65 


2.32 


0.11 


0.04 


-0.0425 


-0.0411 


top 


2.89 


2.66 


2.22 


0.23 


0.005 


-0.0408 


-0.0395 


stacking fftifj^ 


3.20 


2.59 


2.35 


0.13 


0.02 


-0.0417 


-0.0405 


exp (fcc)S'B 


3.27 ± 0.03 2.70 ± 0.03 


2.34 ± 0.02 


0.10 ± 0.03 









TABLE III. Orbital-projected potassium charge for different adsorption sites (5 Ag layers, with substrate rumpling and 
lateral relaxation). 



site 






charge, in |e| 








s 




Py 

{VS X x/3)R30° 


Pz 


total 


fee hollow 


6.513 


4.135 


4.135 


4.052 


18.843 


hep hollow 


6.515 


4.134 


4.134 


4.052 


18.843 


bridge 


6.513 


4.127 


4.144 


4.052 


18.845 


top 


6.502 


4.146 


4.148 


4.051 


18.856 


stacking fault 


6.515 


4.134 


4.134 
(2 X 2) 


4.051 


18.843 


fee hollow 


6.494 


4.115 


4.114 


4.025 


18.757 


hep hollow 


6.496 


4.114 


4.114 


4.025 


18.757 


bridge 


6.496 


4.106 


4.122 


4.027 


18.760 


top 


6.487 


4.123 


4.123 


4.028 


18.772 


stacking fault 


6.496 


4.114 


4.114 


4.025 


18.757 
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TABLE IV. Position of the K 3s and 3p core eigenvalues, and position of the Fermi energy, in (5 Ag layers, with substrate 
rumpling and lateral relaxation). 



site 


(V3 


X V3)R30° 




(2 X 2) 




fee hollow 


-1.253 ; 


-0.659 


-0.069 


-1.256 


-0.662 


-0.062 


hep hollow 


-1.253 ; 


-0.658 


-0.069 


-1.256 


-0.662 


-0.062 


bridge 


-1.253 ; 


-0.659 


-0.069 


-1.257 


-0.662 


-0.062 


top 


-1.253 ; 


-0.659 


-0.069 


-1.258 


-0.664 


-0.062 


stacking fault 


-1.253 ; 


-0.659 


-0.069 


-1.256 


-0.661 


-0.062 



TABLE V. Comparison of the results for the two coverages (fee hollow site). 



structure (V3 X V3)R30°-K (2 x 2)-K 

bond length [A] 3.27 3.20 

effective K radius [A] 1.82 1.76 

binding energy [j^^] 0.042 0.041 

MuUikcn charge on K [|e|] 0.16 0.24 

K 3s, 3p core eigenvalues, relative to Fermi energy [Eh] -1.184 ; -0.590 -1.194 ; -0.600 

work function [En] (clean Ag(lll): 0.131) 0.069 0.062 
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FIG. 1. The structures considered for K, adsorbed on the Ag(lll) surface, at a coverage of one third of a monolayer, 
(\/3 X v^)R30° unit cell. The silver atoms in the top layer are displayed by open circles. The considered poteissium adsorption 
sites arc the top site above the silver atoms with number 1 (filled circles), the threefold hollow sites above atoms 1,2,3 (fee 
or hep hollow, possibly with stacking fault, circles with horizontal lines) or the bridge site above atoms 2 and 3 (circles with 
horizontal and vertical lines). Note that the threefold hollow sites can not be distinguished in this figure. 




X 
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FIG. 2. The structures considered for K, adsorbed on the Ag(lll) surface, at a coverage of one fourth of a monolayer, (2 x 2) 
unit cell. The silver atoms in the top layer are displayed by open circles. The considered potassium adsorption sites are the 
top site above the silver atoms with number 1 (filled circles), the threefold hollow sites above atoms 2,3,4 (fee or hep hollow, 
possibly with stacking fault, circles with horizontal lines) or the bridge site above atoms 3 and 4 (circles with horizontal and 
vertical lines). Note that the threefold hollow sites can not be distinguished in this figure. 
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FIG. 3. Definition of the geometrical parameters. All distances are interlayer distances (definition see text). Note that in all 
the figures, the atoms are drawn purely schematically and the size of the atoms does not scale with the atomic radii. 



d K-Agl 

^1 

dAgl-Ag2 
dAg2-Ag3 =ao/^ 
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FIG. 4. DOS, projected on K, fee adsorption site, (\/3 x y/S)R30° structure. The DOS, projected on all K basis functions, 
is shown, together with the DOS projected on s, px, Py and pz orbitals only. The Fermi energy is indicated with a vertical line. 



Ag (1 1 1) (roots x root3)R30 - K, fee adsorption site 

DOS, projected on K basis functions 




-0.25 



T=0.005 



-0.125 
ENERGY (HARTREE) 
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FIG. 5. DOS, projected on K, fee adsorption site, (2x2) structure. The DOS, projected on all K basis functions, is shown, 
together with the DOS projected on s, p^, Py and pz orbitals only. The Fermi energy is indicated with a vertical line. 



Ag (1 1 1 ) (2x2) - K, fee adsorption site 

DOS, projected on K basis functions 
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FIG. 6. Binding energy -BK/Ag(iii) — ■E'Ag(iii) — Ek, for K adsorbed on the hep site of the (111) surface of a 3 layer Ag slab. 
Various smearing temperatures from 0.001 to 0.05 Eh, and fe-point samplings from 2x2 up to 16x16 have been compared. 



-0.035 



-0.04 



-0.045 
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LU 
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